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We present an investigation of massive star formation that results from the 
gravitational collapse of massive, magnetized molecular cloud cores. We investi- 



gate this by means of highly resolved, numerical simulations of initial magnetized 
. Bonnor-Ebert-Spheres that undergo collapse and cooling. By comparing three 

different cases - an isothermal collapse, a collapse with radiative cooling, and a 
magnetized collapse - we show that massive stars assemble quickly with mass ac- 
cretion rates exceeding 10~ 3 M & yr _1 . We confirm that the mass accretion during 



the collapsing phase is much more efficient than predicted by selfsimilar collapse 



solutions, i.e. M ~ c 3 /G. We find that during protostellar assembly the mass 
accretion reaches 20 — 100c 3 /G. Furthermore, we determined the self-consistent 
structure of bipolar outflows that are produced in our three dimensional magne- 



tized collapse simulations. These outflows produce cavities out of which radiation 
pressure can be released, thereby reducing the limitations on the final mass of 
massive stars formed by gravitational collapse. Moreover, we argue that the ex- 
traction of angular momentum by disk-threaded magnetic fields and/or by the 
appearance of bars with spiral arms significantly enhance the mass accretion rate, 
thereby helping the massive protostar to assemble more quickly. 

Subject headings: accretion, accretion discs, magneto-hydrodynamics, ISM: clouds, 
evolution, methods: numerical 



1. Introduction 



The formation of massive stars is still a highly debated question. Low mass stars accrete 
the bulk of their mass through their circumstellar disks before nuclear burning turns on (eg. 
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Shu et al.l (119871 )) . Massive stars on the other hand, have Kelvin-Helmholtz time scales that 
are smaller than the dynamical time so that young massive stars sta rt to burn their nuclear 
fuels while still accreting the surrounding gas (e.g.. lYorkdl2002l . 120041 . and references herein). 



Infall, or flow of gas through a surrounding disk, therefore faces a major obstacle in the form 
of the radiative pressure that such massive stars will produce as they are still forming. 

Early spherical accretion models suggest that the resulting radiation pressure could limit 
the final mass of the star to ~ 40 if the accretion rate is not high enough ( > 10~ 3 yr" 1 ) 



(l2002h 



( Kahnlll974l ; IWolfire fc Cassinellilll987l ). More recent two dimensional simulations by lYorke fc Sonnhalter 



showed that the limitations on the final mass of the massive star can be relaxed if 
accretion through the protostellar disk is inclu ded in models of massi ve star formation. Still, 
even this result yields mass limits of ~ 43 M . iKrumholz et al.l (120051 ) argued that the effects 



of radiation pressure are limited by the escape of radiation through an outflow cavity. These 
results are were based on simulations using Monte-Carlo-diffusion radiative transfer models 
where the outflow cavities are parameterized by varying opening angles. 



Alternatively, massive stars could form through coalescence of intermediate mass stars (IBonnell et al 



19981 ). This formation process would be starkly different from the formation of low mass 
stars which assemble quickly through accretion of the molecular gas. So far, observations of 
the intermediate state of massive star formation are rare and difficult to obtain. Neverthe- 
less, a fe w massive objects which show evidence for an ongoing accretion process a re known 
by now (jchini et alibooi IPatel et al.lhood : Ichini et al.lbood : beltran et al.ll20Q6h . 



We have developed detailed 3D collapse simulations of magnetized, star forming cores 
and have shown that outflows and gravitational collapse are inseparably linked. In fac t, we 
found in earlier work on the assembly of low mass stars (e.g.. iBanerjee fc Pudritzll2006l ) that 
outflow cavities are carved out of the collapsing envelope very early - much before even a solar 
mass of material has been assembled into the star. This important early outflow activity we 
suggested, could play a key role in actually opening up a channel to radiation pressure in 
the more massive systems. The outflow channel is not imposed by an external model but 
arises as a natural consequence of the generation of an outflow during the formation of a 
magnetized disk. 

Our first study of massive star formation was in the context of the turbulent fragmenta- 
tion picture. Specifically, we followed the formation of massive stars in supersonic turbulence 
by using a 3D, adaptive mesh (AMR) simulation and showed that filamentary accretion 
plays an important role in foc using the collapse of material onto a massive young star and 
its disk (IBanerjee et al.ll2006l . henceforth BPA06). We included a full list of coolants and 
found that the central regions of the collapse behaved somewhat like Bonnor-Ebert spheres. 
We did not have MHD at work in those calculations. 
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The present paper follows up on these previous efforts. Instead of trying to track the 
formation of a massive star in a very large cluster simulation, we concentrate here on the 
collapse of one massive central region chosen to mimic the massive collapse we found in 
BPA06. The choice of an initial magnetized Bonnor-Ebert sphere is a reasonable starting 
point for such a focused-in study. In particular, we study the collapse of a massive molecular 
cloud core with a total mass of ~ 170 M , which we model as a sphere that is initially in in 
pressure equilibrium with the surrounding medium i.e. Bonnor-Ebert-Spheres (BE-Spheres). 

We shall show that outflows are also an early aspect of the formation of massive stars. 
We self-consistently determine the initial scale and structure of outflow cavities that are 
carved by magnetic tower flows that are driven by the forming protostellar disk. In particular 
we study the influence of radiative cooling and magnetic fields on the collapse of such cloud 
cores thereby discussing the similarities and differences of these models. We find that these 
outflows are launched very early, much before the central stars have even accumulated a 
solar mass worth of material. This means that even before radiation pressure becomes a 
factor in the formation of a massive star, the outflow has created the cavity out which the 
radiation pressure, when it eventually appears, can be released. Early outflows therefore, 
play a central role in the assembly of such stars. 

Our paper is organized as follows. In Sec. [2] we describe the numerical modelling and 
initial conditions of our simulations. We discuss the accretion process and rates in Sec. [3j 
In Sec. H] we discuss the angular momentum evolution and distribution for the different 
models, and we discuss the influence of the early outflows in Sec. [5J Finally, we conclude the 
results of this work in Sec. [6j 



2. Numerical modeling 



For this study we follow the collapse of Bonnor-Ebert- type (Ebert 



molecular cloud cores using the AMR based FLASH code (IFryxell et al 



1955 



Bonnorlll956f ) 



2000) to solve- tho 



gravito-magneto-hydro dynamic equations (de tails on the Bonnor-Ebert-Profile and of this 
setup can be found in iBanerjee et al.l (12004 ). henceforth BPH04). We study the collapse 
of three different cases which all have the same initial density and temperature distribution. 
The three cases are: 



The isothermal case (iso): Here, we assume an isothermal equation of state (adiabatic 
index, 7 = 1.0001) throughout the calculation. 

The hydro case (hydro): For this case we incorporate the cooling ability from molecular 
and dust radiation in the optical thin regime and use a radiation diffusion approxima- 
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tion in the optical thick regime. Details of our cooling treatment can be found in the 
appendix of BPA06. 

• The magnetized case (mag): Here, we apply an initial magnetic field parallel to the 
rotation axis (z-axis) and use the same cooling as in the pure hydro case. 



The three different cases presented in this study, have the same initial values for the total 
mass of the cloud core, M = 168 M , its radius, R = 3.34 x 10 5 AU = 5 x 10 18 cm, and the 
BE-core density, po = 3.35 x 10 -21 gcm -3 . We also enhance the density of all cores by 10% 
and add a 10% m = 2 density perturbation. The initial free fall time associated with the core 
region is tg = a/37t/32 G po = 3.63 x 10 13 sec = 1.1 x 10 6 yrs. We also give the cloud cores 
a slight spin with a constant angular velocity. The two non-magnetic simulations start with 
an angular velocity of Q — 5.5 x 10~ 15 rad sec -1 whereas the magnetized case rotates with 
an angular velocity of Q — 1.1 x 10~ 14 rad sec -1 . We increase the initial rotation by a factor 
of two in the magnetize d case because magnetic b raking spins down the sphere substantially 



prior to the collapse (see lBanerjee &: Pudritzll2006l ) . These initial conditions are equivalent to 



an fltg of 0.2 and 0.4 for the non-magnetized and magnetized cases, respectively. In BPH04 
we showed that the former case results in a ring which fragments into a binary system if 
only molecular cooling is considered. The initial isothermal sound speed, Ci SO , is in all cases 
0.408 km sec -1 . 

Initially, the magnetic field in the magnetized case is given only a component parallel 
to the rotation axis (z-axis) where B z varies between 0.33 /iG and 1.36 fiG corresponding 
to a constant (3 = p/(B 2 /8ir) = 76 in the equatorial plane. This gives a mass-to-flux 
ratio of M/$b = 3.2/yG which is highly supercritial (the critial value is ~ 0.12/ y/G). 
Nevertheless the magnetic field becomes dynamically important during the collapse and will 
be able to drive outflows and jets as we will see later . This choice of a supercritical core 



is in agreement with recent work by iTilley fc Pudritzl (120051 ) which studied the formation 



of turbulent fragmentation and the IMF in magnetized turbulent clumps. In that paper, it 
was demonstrated that substantially supercritical clumps did a far better job in reproducing 
the observations of the core mass function in cluster forming environments, than did clumps 
that were nearer to magnetically critical values. One way that this can easily occur in a 
turbulent environment is for shocks to rapidly assemble cores by flows along field lines. 

The choice of our initial c ore mass is taken from recent observations of cores in regions of 



high mass star formation (e.g. jReid &: Wilsonll2006l ). The actual physical values we use arise 
in this context arise from the consideration of a BE sphere with a core density of 10 3 cm -3 
(3.35 x 10 -21 gcm~ 3 ) at 20 K which is a typical value of such molecular cores. As we will 
see later, the mass accretion rate depends only weakly on the initial mass of the core. 
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3. High accretion rates during the collapse phase 



The collapse of spherical cloud cores 



alytically and numerically (see e.g.. lLarson 



Whitworth fc Summers 1985; Foster & Chevalier 



las been studied by many aut 



1969 



Penston 



1993 



1969 



Shu 



Hennebelle et al. 



l ors, bot 



1977 



2003 



Hunter 



h an- 



1977 



Banerjee et al. 



20041 ). In the case of an initial singular isothermal sphere (SIS) there exists an elegant self- 
similar solution to this problem wherein an expansion wave travels from the (singular ) center 
outward with the speed of sound, thereby initia ting an inside-out collapse (jShulll977l ). How- 
ever, as pointed out by IWhitworth et al.l (119961 ) a SIS configuration is unnatural bec ause - 
among other difficulties - its collapse can not produce binaries (see also |Pringldll989l ). The 
collapse of non-singular (ie Bonner -Ebert ) cores proceed differently than singular spheres 
and have distinguishable implications. Firstly - as demonstrated by many authors (e.g., 



Larsonlll969l ; |Penstonlll969l ; iFoster &: Chevalierlll993l ; iHennebelle et al.ll2003l ; iBanerjee et al. 



2004 ) - the collapse proceeds from the outside-in rather than from the inside-out and the 
density maintains a flat profile at the core center, where the core size is of the order of the 
local Jeans length at every epoch (see also IWhitworth fc Summers! Il985l . for a summary of 
analytic solutions). Secondly, the radial distribution of the infall velocity peaks at the edge 
of the flat density core and falls off quickly towards the center (e.g. see Fig. 3 in BPH04). 
The velocity also becomes supersonic as the core density increases and the size of the flat 
region shrinks. Recent observations of a pre-Class object show that the collapse proceeds 
supersonically, stro ngly supporting a La rson- Penston-type collapse rather than an expansion 
wave-type collapse (IFuruya et al.ll2006l ). Thirdly - and of great importance - the mass ac- 
cretion in the non-selfsimilar case is much higher than predicted from the selfsimilar collapse 
of a singular isothermal sphere. We find that the mass accretion in the early phase of the 
collapse is 

c 3 

Afw 20 - 100— (1) 

Gr 



(c is the isothermal sound speed and G is Newton's gravitational constant). Note the self- 
similar SIS collapse gives a mass accretion of only 0.96 c 3 /G. Typical values of the sound 
speed in cold cloud cores (T ~ 20 K) are of the order of a few 10 2 m sec -1 which gives 
c 3 /G ~ 10~ 6 — lO~ 5 M yr -1 . Our res u lt fro m n umerical s imula tions are in agreement 
with the early analytic result s of lLarson fjl969h and iPenstonl (119691 ) (see also iHunterl 119771 ; 
Whitworth fc Summers! Il985l ). 



The remarkable point of this result is that the high accretion rates are achieved even 
without initial turbulence and during the isothermal phase of the collapse. The main reason 
for the high accretion rate in this idealized case is the supersonic infall velocity, v r , close to 
the peak density. Even moderate Mach numbers, A4, of 2 — 3 (The Larson asymptotic Mach 
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Fig. 1. — Shows the time evolution of the column density and mass accretion in the mag- 
netized case. We assume the beginning of the collapsing phase at t = to when the column 
density reaches S core « lgcm -2 and the pressure is P/ks ~ 10 8 Kcm~ 3 (see Fig. [2]). At 
t = to the mass density in the core is p core = 2.44 x 10~ 17 gcm -3 corresponding to a free fall 
time of tfi = 1.35 x 10 4 yrs. At t — t + 1.1 tg the core density reached already 10 4 gcm~ 2 . 
The initial profile of our simulation is marked by t = 0. 




1 = 

t = t + 14201 y 
t = t + 14714 y 
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Fig. 2. — Shows the time evolution of the pressure and temperature in the magnetized case. 
See Fig. [1] for the corresponding column densities and label explanations. The pressure scales 
as P oc E 2 and the temperature develops a R~ os profile within the warm, dense region. 



- 7- 



number is M. ~ 3) enhances the mass accretion rates relative to c 3 /G because 

c 3 /G^v*/G = M 3 c 3 /G 



(2) 



in the supersonic limit. Additionally, the core is continually embedded in a high pressure 
environment. 

One-dimensional numerical studies o f the long term evol u tion ( i.e. beyond one dynam- 
ical time) of a collapsing BE- Sphere by iFoster fc Chevalier! (119931 ) showed that the mass 
accretion is not constant but decreases with time after reaching a peak value. Beyond this 
peak, the remaining gas in the envelope drains out with an ever smaller accretion rate over 
time. This situation could be different for the collapse of realistic cloud cores, which are 
certainly not isolated objects but are instead surrounded by a clumpy media and attached 
to filaments wherein accretion might sustain high rates for a longer time. 

Using the magnetized simulation as the fiducial case, we show the evolution of the col- 
umn density and corresponding mass accretion rates in Fig. [T] and in Fig. [2] the evolution of 
the pressure and temperature as functions of the disk radius 0. When the density approaches 
£ ~ lgcm -2 (we denote the corresponding time as to) the mass accretion rate is already 
~ lO~ 4 M yr~ 1 and the surrounding p r essur e is 10 8 K cm" 3 ks- These results are in agree- 



ment with the model of iMcKee fc Tan! (120021 . 120031 ) who showed that this high pressurized, 
compact cloud cores accrete gas with accretion rates as high as 10 _3 MQyr _1 . We find that 
accretion rates of this order (~ 10~ 3 M & yr _1 ) are reached within only ~ 14.7 x 10 3 yrs which 
corresponds to 1.1 dynamical times. The core column density after one dynamical time is 
~ 10 4 gcm -2 while the pressure achieves a value P/kg « 3.5 x 10 16 Kcm~ 3 . As long as the 
core remains isothermal (the most efficient cooling regime), the pressure scales with the col- 
umn density as P oc £ 2 (i.e. the pressure profile is close to R~ 2A ) and slightly steeper in the 
inefficient cooling regime where the temperature rises during the collapse (see right panel of 
Fig. [2]). The continuously high external pressure and the supersonic infall velocity maintains 
the high mass accretion. We also show the mass accretion scaled to the quantity cf so /G (c iso 
is the initial isothermal sound speed) in Fig. [TJ From this one can see that accretion is much 
more efficient (by a factor of ~ 100) than expected from the collapse of a SIS. 

These results show that high mass accretion rates are a natural result of the early collapse 
phase of (non-singular) collapsing cloud cores with flat-topped density profiles. Turbulent 
driving for the rapid assembly of massive stars is not a necessary ingredient, although it doe s 
further enhance the accretion rate as shown in analytic models ( IMcKee &: Tan! 120021 . 120031 ). 



1 All quantities which are functions of the disk radius, R, are azimuthally averaged and density weighted, 
e.g. f(R) = (2nY,(R)y 1 J dzd<pp(x) /(x), where E(i?) = {2-k)- 1 J d<pdzp(x) is the column density. 
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The enhancement is even larger in fully realized 3D turbulence, due to filamentary accretion 
(see e.g., BPA06). 

For comparison, we summarize the collapse phase of the three different cases (mag, 
hydro, iso) in the Figs. [3] -[5] at the time when the column density reaches 5 x 10 3 gcm -2 . 
Common features can be seen mainly in the envelope of the collapsing cores. All three cases 
end up with a column density profile which is close to oc R~ 12 in the disk envelope (see left 
panel of Fig. [3]). In the inner part of the disk (< 100 AU) the influence of the (inefficient) 
cooling processes and magnetic fields become visible. The non-isothermal cases develop a 
steeper density profile as infalling material tends to pile up at the shock regions at the edge 
of the warm core. These shock regions appear at a radial distance of ~ 20 AU because the 
core becomes opaque inside this radius (BPA06). The difference with the isothermal case 
inside < 100 AU can also be seen in the mass distribution of the disk which we show in the 
right panel of Fig. [3j More material accumulates in the hydro and magnetized cases at the 
time when all three cases reached a peak column density of 5 x 10 3 gcm -2 . 

We show the comparison of the mass accretion rates of the three different cases in Fig. |H 
The right panel gives the mass accretion in terms of solar masses per year and compared to 
the initial isothermal sound speed, whereas the right panel shows this value scaled to the 
local sound speed quantity cf ocal /G. In physical terms (i.e. M /yr) the hydro case is most 
advanced and accretes gas with several lO _3 M yr _1 corresponding to ~ 300cf so /G. Then 
again, the accretion rates of all three cases are very similar when compared to local quantity 
c 3 ocal /G (~ 40 cf ocal /G). This shows, again, that the estimate for the mass accretion in the 
supersonic regime is sensitive to the Mach number and the local sound speed. To complete 
this picture, we show the infall velocity in Fig. [5] measured in km/ sec and Ci SO (left panel) 
and in terms of the local sound speed (right panel). The high infall velocity (and therefore 
the high mass accretion) in the hydro case is due to the warm core region which has a high 
sound speed. The local Mach number of all three cases differs only within a range 2 — 2.5 
whereas the isothermal Mach number ranges from 2.5 — 4 between them. 

We also point out that the mass accretion during this early stage of collapse does not 
strongly depend on the initial mass of the c loud core (the good estim ate of Eq. [2] does 



not have a mass dependence). In BPH04 and lBanerjee fc Pudrita (120061 ) we showed that a 



2M cloud core reaches also mass accretion rates of 10 4 — 10 3 M Q /yr. 
In summary, we find that the typical infall velocity of a collapsing cloud co re is 2 — 3 q 



which recovers earlier re sults studying the isothermal collapse phase (e.g., iLarsonl 11969 



Foster &: Chevalieri Il993l ). The supersonic infall velocity is the main reason for the high 
mass accretion of collapsing cloud cores. Furthermore, the increased local sound speed due 
to thermal effects in the warm core region enhances the mass accretion in the post-isothermal 
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Fig. 3. — Comparison of the column density profiles and disk masses for three different cases, 
hydromagnetic (mag), pure hydro (hydro), and pure isothermal (iso) simulation at the time 
when they reached the same core column density S core ~ 5 x 10 3 gcm~ 2 . 
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Fig. 4. — Comparison of the mass accretion rate for the three different cases, hydromagnetic 
(mag), pure hydro (hydro), and pure isothermal (iso) simulation at the time they reached the 
same core density (see Fig. [3]). The left panel shows the mass accretion in units of M yr _1 
whereas the right panel shows the same quantity compared to c 3 ocal /G, where ci oca i is local 
sound speed. 
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regime. 



4. Angular momentum extraction by bars and magnetic fields 

Most stars, if not all, rotate. The origin of stellar angular momentum is in itself an 
important problem in star formation theory. It has a natural explanation in the turbulent 
fragmentation picture, where it arises from the oblique shocks that create the dense cores 
(eg. Tilley & Pudritz 2004). Strict conservation of the specific angular momentum in a 
core would inhibits accretion through the protostellar disk whenever the infalling gas hits its 
centrifugal barrier. The issue of redistributing or extracting the angular momentum within 
the protostellar disk is important for understanding the assembly of (massive) stars. 

In Fig. [6] we show the evolution of the specific angular momentum in the magnetized case 
as a function of radius (left panel) and enclosed mass in the disk (right panel). Infalling, high 
angular momentum gas increases the specific angular momentum at small radial distances 
with time. Apart from the advection of angular momenta due to the accreting gas, there is 
additionally, a resorting of disk angular momentum operating due to turbulence. This can 
be seen from the right panel of Fig. [6] where we plot the specific angular momentum as a 
function of the enclosed mass. This is a useful thing to do because this quantity is constant 



with time if angular momentum is not redistributed within the disk (see e.g., lAbel et al. 



20021 ). For the magnetic case shown in Fig. [6] the angular momentum is decreasing with time 
within the disk (inside Mdi S k|$ 1M S ). This decrease of angular momentum at a given mass 
shell allows for continuing efficient accretion through the disk plane. In particular, we find 
that the specific angular momentum changes by an amount of Aj z ~ 5.3 x 10 19 cm 2 sec -1 (i.e. 
by a factor of 2.8) at a mass shell of 0.1 M & within 14.7 x 10 3 yrs. In the magnetized case the 
angular momentum is extracted by magnetic braking, magnetic torque, the launching of an 
outflow, and by the development of a bar in the disk. Depending on the field geometry and 
strength magnetic braking b y a poloidally dominated magnetic field can be an efficient way to 



extract angular momentum (IMouschovias &: Paleologoulll979l . 119801 ) . An additional toroida l 



component of the magnetic field gives rise to a magnetic torque (e.g.. lKonigl fc Pudritzll2000l ) 



which in the case of rotationally generate d toroidal field (backwards bend field) extracts 



angular momentum from the protodisk. In iBanerjee fc Pudritzl (120061 ) we showed that the 



magnetic torque can be a substantial fraction (10% - 50%) of the angular momentum flux 
from the accreting gas. 

In Fig. [7| we compare the angular momentum distribution and the rotational velocity 
for the three different cases (Note that the overall enhancement of the angular momentum 
and rotational velocity in the magnetized case is due to the initially higher angular velocity). 
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Fig. 5. — Comparison of the radial infall velocity for the three different cases, hydromagnetic 
(mag), pure hydro (hydro), and pure isothermal (iso) simulation at the time they reached 
the same core density (see Fig. [3]). The left panel shows the infall velocity in km/sec and 
compared to the initial isothermal sound speed, Ci SO , and the right panel shows the same 
quantity measured with the local sound speed, ci oca i - 
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Fig. 6. — Shows the time evolution specific angular momentum as a function of radius and 
enclosed mass in the magnetized case. See Fig. [1] for the corresponding column densities and 
label explanations. 
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Fig. 7. — Comparison of the specific angular momenta and rotational velocity for the three 
different cases, hydromagnetic (mag), pure hydro (hydro), and pure isothermal (iso) simula- 
tion at the time they reached the same core density of 5 x 10 3 gcm -2 (see Fig. [3]). The left 
panels show this quantities as a function of the radial distance and the right panels show the 
same quantities as functions of the enclosed disk mass, M enc i. Note that the initial rotation 
in the magnetic case is two time larger than in the two other cases (see also Sec. |2]) 
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During this early stage of the collapse, the disk is still more massive than the central object (s) 
and is still out of equilibrium. In all three cases the rotational velocity is sub-Keplerian with 
10% to 60% of the centrifugal velocity Apart from the general common features of the 
angular momentum distribution of the three models there are important differences. Both of 
the non-isothermal cases show a 'pile-up' of high angular momentum gas at the shock fronts 
which sets the edge of the disk. This can be seen in the lower left panel of Fig. [7| where 
the high rotational velocity material accumulates at the shock distance of few x 10 AU. 
Inside the protodisk (inside the shock fronts) the angular momentum is very efficiently 
transported by the large bars and/or the magnetic field (see Figs. [8] and [9]). Therefore, 
the rotational velocity drops rapidly towards smaller radii in the hydro and magnetized case. 
The isothermal case shows a much more evenly distributed toroidal velocity when compared 
to the Keplerian velocity. The reason here is the absence of shocks in the disk plane which 
allows a continuous inflow of high angular momentum gas. Additionally, the absence of a 
strong bar in the isothermal case does not allow for efficient transport of angular momentum 
(see Fig. CUD. 

We investigate this further by presenting 2D snapshots of the disk planes for the three 
different cases in Figs. E] - DZ3 on a scale of ~ 300 AU and ~ 40 AU, respectively. Both the 
magnetized (Fig. [8]) and the hydro simulation (Fig. [9]) developed a large bar with an extent 
of ~ 250 AU. Comparing the small scale structure of the magnetized and hydro case (right 
panel of Figs. El and [9]) one sees that the former case develops a larger, more rotationally 
supported protodisk. This is due to the additional magnetic pressure and the higher toroidal 
velocity in the magnetized case. Both cases accrete efficiently gas through their bars/spiral 
arms onto the protodisk. This can be seen particularly well in the right panel of Fig. [9] 
where high density gas is streaming along the spiral arms towards the central object. The 
isothermal case also developed a bar but it is much smaller and less pronounced. Therefore, 
the drop of the rotational velocity (measured in Keplerian velocity) towards the disk center 
is only a factor of two (in the hydro case the rotational velocity decreases by a factor of 6). 



5. Early outflows during massive star formation 

Little is known observationally about the influence of outflows on the early assembly 
of massive stars. On one hand, they could reduce the mass accretion onto the massive 
(proto)star if the outflow carries a substantial mass. On the other, early outflows provide a 
natural anisotropy of the accreting gas which results into low density cavities. Such cavities 
are like funnels out of which the radiation from the already active star can escape. Without 



such a pressure-release valve, trapped radiative flux would halt the in-fall. iKrumholz et al. 
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Fig. 8. — Magnetized case: The panels show 2D slices through the disk plane at two different 
scales (large scales at right and small scales at the left). This snapshot is taken 1.4 x 10 4 yrs 
into the collapse, corresponding to S coro « 5 x 10 3 gcm~ 2 (see Fig. [3]). 
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Fig. 9. — The hydro case: the panels show 2D slices in the disk plane at two different scales 
(large scales at right and small scales at the left). This snapshot is taken 1.4 x 10 4 yrs into 
the collapse, corresponding to S core ~5x 10 3 gcm~ 2 (see Fig. EJ). 
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Fig. 10. — The isothermal case: the panels show 2D slices in the disk plane at two different 
scales (large scales at right and small scales at the left). This snapshot is taken 10 4 yrs into 
the collapse, corresponding to S corc « 5 x 10 3 gcm~ 2 (see Fig. [3]). 




Fig. 11. — Shows 2D snapshots perpendicular to the disk plane of the magnetized case 
at different times on large scales (~ 300 AU). The left panel shows the situation at t = 
to + 1.31 x 10 4 yrs (0.98 iff) into the collapse and the right panel shows the configuration 
1566 yrs later after magnetic pressure inflates the disk and the outflow is launched (see also 
Fig. d2]for a close-up). 
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( 120051 ) studied the effect of outflow-funnels using a Monte Carlo radiative-transfer method 
which shows that radiation pressure is greatly reduced by radiation escaping through the 
outflow cavities. 

Magnetic fields coupled to the protostellar disk can be the driving power for such out- 
flows. A variety of self-consistent simulations of collapsing magnetized cloud cores show 
that outflows are launched if t he toroidal magnetic field pressure overcomes the gravi 



Banerjee fc Pudritz 



tational force (Tomisakal Il998l . 12002 ; Matsumoto fc Tomisakal 12004 : Machida et all 12004 



2006) . This early- type outflows can be understood in terms of a grow- 



ing magnetic tower (jLynden-Bellll2003l ) : The rotating (proto)disk generates a strong toroidal 



field component by winding up the threading field lines. The resulting magnetic pressure is 
in local equilibrium with the gravitational force and the ram pressure of the infalling mate- 
rial. But every further rotation increases the toroidal field component thereby shifting the 
equilibrium location (which is characterized by a shock) towards higher latitudes. The result 
is an inflating magnetic bubble in which material is lifted off the disk. 

Some of the highlights of this paper are shown in the following figures, wherein we 
explicitly demonstrate the launch of the outflow. We present 2D cuts through our 3D data, 
which are shown as the images in Figs. [TTJ and [12j These show the collapse state at different 
times and physical scales. The time sequence shown in Fig. [TTJ demonstrates the situation 
shortly before the launching of the outflow (left panel) and roughly 1600 yrs later (right 
panel). On this large scales (~ 300 AU) one can see the strong bending of the dragged-in 
magnetic field lines (the toroidal component is not shown in the 2D slices) and the inflated 
protodisk. Also clearly visible are the shock fronts which separate the disk from the accreting 
environment. 



Close-ups (~ 20 AU) of the early stage of the outflow launching are shown in Fig. 
where we zoom in by a factor of 10 in comparison with the previous figures. Here again the 
collapsing stage shortly before the flow reversal is shown in the left panel and the onset of the 
outflow (less than 200 yrs later) is shown in the right panel. Due to the stronger magnetic 
field strength deeper in the gravitational potential well, the outflow velocity (~ 4 km sec -1 ) 
is higher than on larger scales which powers a stronger outflow. 

It is also worth emphasizing that outflow only occurs when magnetic fields are included 
in the simulations. In this paper (as in our others), we have never seen outflows associated 
with our purely hydrodynamic collapses. We conclude that purely hydrodynamic collapse 
misses a critical ingredient in stellar formation - namely- the early launch of disk related 
outflows. 
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6. Summary and Conclusions 



In this work we simulated the 3D collapse of magnetized, massive molecular cloud cores 
which are the progenitors of massive stars. We carried out this investigation at early stages 
of the collapse, long before a massive star has even assembled at the center of the resulting 
disk. Using three different 3D collapse simulations (an isothermal case, a pure hydro case 
with radiative cooling, and a magnetized case) of Bonnor-Ebert type cloud cores we show 
that the mass accretion during the collapse phase is much more efficient than calculated 
from the collapse of a singular isothermal sphere (SIS). 

We find that the mass accretion rates can reach ~ 10 2 c 3 /G which is a factor of 10 2 higher 
than the SIS case. This high mass accretion is due to the supersonic infall velocity and the 
high pressure at the edge of the collapsing core, where the pressure increases with density as 
P oc X 2 . The mass accretion can be estimated by M ~ v f/G = M 3 c 3 /G where the Mach 
number, A4, can reach 3 or more with infall velocities ranging from 1 — 1.5 km sec -1 . Note 
that, due to inefficient cooling at densities of S ^ 10 g cm -2 the temperature and therefore the 
sound speed increases which in turn results in a higher physical mass accretion rates. There- 
fore, mass accretion rates of 10 -3 M yr" 1 and higher are easily reached in the early stage of 
the assembly of massive stars. This high accretion rates allow the formation of massive stars 
in only 10 4 yrs. A formation time scale of only 10 4 yrs is short enough for the massive star 
to build up before its surrounding disk is p hoto-evaporated from the UV environment emit - 
ted from the surrounding OB association (iHollenbach et al.ll2000i ; iRichling fc Yorkdll998l ). 
This allows also for an unified picture of star formation where low and high mass star form 
the same way from the collapse of an overdense cloud core. This unified theory is strongly 
supported by recent obs ervations of the assembly of massive stars through accretion disks 



Jchini et allbooi hood) a nd fro m the detection of high accretion rates onto a very young 



massive star (IBeltran et al.ll2006l ). 



Furthermore, we find that angular momentum is efficiently transfered by bars/spiral 
arm and magnetic torque. The simulations with radiative processes develop a larger bar-like 
structure than the isothermal case and are less rotationally supported in the inner region 
of the disk. The reduced specific angular momenta also enhance the accretion rate as low 
angular momentum gas can be accreted more efficiently In particular we showed that the 
specific angular momentum decreases by a factor of 2.8 within 14.7 x 10 3 yrs (1.1 tg) into the 
collapse. In all cases, an equilibrium disk did not yet develop because at this early stage, 
the disk is still much more massive than the central object. At this stage the rotational 
velocities are sub-Keplerian and and range from 0.4 — 0.7t>K C p within the protodisk. 

Outflows play an important role in this process by reducing the radiation pressure which 
can escape through funnels carved by the expelled gas (IKrumholz et al.ll2005l ). Here, we have 
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shown that that outflows are a natural consequence of the early stages of collapsing, massive 
magnetized (and not purely hydrodynamic) cloud cores. This process is not restricted just to 
the case of low mass star formatio n which we prev iously examined. Recent observations indi- 



cate also early type outflows (e.g.. IShepherdll2003l ) which might be associated with magnetic 
fields. 

The critical issue that we have uncovered is that for both low (our previous work) and 
high mass star formation (this work), magnetic pressure from a wound-up toroidal field 
becomes eventually strong enough to reverse the infalling gas and launch an outflow. We 
find that the outflow has the shape of a hollow cylinder which o utflow speed reaches up to 



5 km sec 1 at 20 AU above and below the disk midplane (see e.g. iBanerjee fc Pudrita 12006 
on details on magnetically driven outflows). 

We conclude that this universality in the formation of disks and the early formation of 
outflow cavities that are created by the jets that are launched from their surfaces, implies 
that the radiation pressure from massive stars appears too late on the scene to substantially 
affect the assembly of a massive star by disk accretion. We are pursuing more detailed 
computations of this picture that will eventually include radiative feedback within more 
complex, cluster forming environments. 



We thank the anonymous referee for some useful remarks. The FLASH code was de- 
veloped in part by the DOE-supported Alliances Center for Astrophysical Thermonuclear 
Flashes (ASC) at the University of Chicago. Our simulations were carried out on computer 
clusters of the SHARCNET HPC Consortium. REP is supported by the Natural Sciences 
and Engineering Research Council of Canada. 
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Fig. 12. — Close-up 2D snapshots of the disk region (perpendicular to the disk in the yz- 
plane) in the magnetized case. The left panel shows the situation at t — 1.45 x 10 4 yrs (1.08 tg) 
into the collapse and before the flow reversal and right panel shows the configuration 188 yrs 
later when the outflow is clearly visible. 
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